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Abstract 

We shall construct a class of nonlinear reaction-diffusion equations 
starting from an infinitesimal algebraic skeleton. Our aim is to explore 
the possibility of an algebraic foundation of integrability properties and of 
stability of equilibrium states associated with nonlinear models describing 
patterns formation. 
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1 Introduction 

In his famous paper, Turing suggested that a system of chemical substances 
reacting together and diffusing through a tissue, could describe the main phe¬ 
nomena of morphogenesis m ■ In particular, in his work it was emphasized 
that patterns could appear if one of the substances diffuses much faster than 
the other. 

Nonlinear reaction-diffusion systems have then been proposed to answer the 
question about how cells, under the influence of their common genes, could 
produce spatial patterns, see e.g. [J| and references therein. They consist in 
models describing generation of patterns from an initially homogeneous state 
taking into account the relevance of chemical gradient in biological systems, in 
particular nonlinear interaction of two chemicals and their diffusion. Patterns 
formation turns out then to be the output of local self-enhancement such as 
local autocatalysis and long range inibition. A simple model proposed by Koch 
& Meinhardt is the activator-substrate systems (constituted by a self-enhanced 
reactant and a depleted reactant which plays the role of the antagonist). In 
the simplest mathematical form, only few relevant parameters characterize the 
model: the normalized diffusion constant and the normalized cross-reaction 
coefficient. In particular, the inibition due to the substrate reactant can be 
effective if the normalized diffusion constant is << 1 (that means the diffusion 
constant of the activator should be much lesser than the diffusion constant of 
the substrate), a necessary condition for the generation of stable patterns. 
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Recent studies pointed out how, beside being a universal principle explain¬ 
ing regular pattern formation in chemical, physical and biological morphogenesis 
systems, Turing’s activator-inibitor principle is also at the basis of regular pat¬ 
tern formation in a variety of ecosystems. In particular, the scale-depending 
feedback is a unifying ecological principle playing a key role and mainly con¬ 
sisting in short-range ecological facilitation such as local modification of the 
environment and long-range competition for resources; see e.g. na whereby the 
prerequisite of long-distance negative feedback is proposed as a basic principle 
for regular pattern formation in ecosystems. 

In this paper we show that this condition can be recognized as a condition on 
the internal symmetry properties of a system. This perspective, concerned with 
the action of a group of transformations on the space of possible configurations 
(‘fields’), will be exploited in Section [2] and detailed in the Appendices. 

In particular, in Section [3] we shed new light on the long-range negative 
feedback condition by a case study exploiting the internal symmetry algebra 
of twisted reaction-diffusion equations with a null basic production term. We 
observe that, in the limit of a null normalized diffusion constant, such a condition 
is related with the existence of soliton solutions - i. e. travelling waves remaining 
stable after interactions. We can see that requiring the normalized diffusion 
constant being null corresponds to look for specific symmetry properties of the 
associated system; indeed this integrable case can be seen as a limiting case of 
the internal symmetry algebra (see Appendix 2). 

As a major result, in subsection 13.21 we recover activator-substrate systems 
by performing a slight modification of the internal symmetry algebra of twisted 
reaction-diffusion equations. We obtain as a byproduct that the necessary condi¬ 
tion for the generation of stable patterns for such a system, besides being related 
with general integrability properties in the limit of a null normalized diffusion 
constant, can also be formulated in terms of ‘closeness’ properties within the 
symmetry algebra vector space. 

This moreover suggests the possibility to enhance the interplay of the two 
different approaches to the study of stability: the qualitative study based on 
the analysis of the stability of equilibrium configurations and global integrability 
properties related with existence of stable solutions such as travelling waves or 
solitons. 

The scope of this work is therefore to propose an intrinsic algebraic form of 
the local mechanisms expressed by reaction-diffusion partial differential equa¬ 
tions. Our considerations are based on the well known duality between lin¬ 
ear differential forms and tangent vector fields on manifolds, and make use 
of algebraic and geometric techniques developed within a theoretical physics 
framework, see e.g. n nu in in m, in particular suitable generalization of 
the structure equations of a Lie group. In the present note, we show that the 
appearing and the significance of parameters in a model can be characterized 
within an algebraic-geometric formulation in terms of integrable “towers with 
infinitesimal algebraic skeletons” (see Appendix 1), which generalizes that of 
integrable connections on spaces of configurations with symmetries. 
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2 Underlying algebraic skeletons 

Concerning real ecosystems in homogeneous landscapes, establishment and sur¬ 
vival of organisms can be inhibited by limited resources or other stress factors: 
it is an empirical evidence that, although the organisms can acts only locally, 
the effects have influence at distance and therefore are of global nature, showing 
a dependence on density of the organisms themselves Ca¬ 
in this paper we propose an approach which is based on the study of global 
properties of partial differential equations such as internal symmetries and in¬ 
variance properties having however an issue in dynamics. The underlying idea 
is that transformations of configurations of a system can be globally studied 
by means of the theory of the action of Lie groups on manifolds. Differential 
equations therefore are an issue of the differential content carried by a Lie group 
(and its Lie algebra) and by its structure equations providing connections on 
the space of configurations. 

Although such mathematical tools could appear of a quite abstract nature, 
nevertheless they are the right tools to deal with global properties at large 
scales, being natural tools for connecting local data to global ones. In particular 
they implement the concept of ‘positional information’ m- In Appendix 1 the 
general context is sketched and appropriate references to the abstract theory 
are given. 

It is well known that, looking for invariance properties, a “prolongation” 
algebra can be associated with a given system of nonlinear partial differential 
equations. Nonlinearity, then, results in a quite intriguing algebraic structure 
which is “only partially” a Lie algebra (we speak of an ‘open’ Lie algebra struc¬ 
ture) [12] • By an inverse procedure based on the intrinsic duality between Lie 
algebras and differential systems [ 2 ], open Lie algebraic structures can ‘gener¬ 
ate’ whole families of different nonlinear systems bound by the same internal 
symmetry structure. 

We show that a slight modification of the internal symmetry properties gen¬ 
erates new models which can contain possible integrable subcases. As an ex¬ 
ample, by a modification of the symmetry algebra associated with a model for 
pattern formation on the shells of molluscs [8] we shall recover the activator- 
substrate reaction-diffusion model proposed by Koch & Meinhardt. In both 
models the necessary condition for the generation of stable patterns, i.e. the 
prerequisite of long-distance negative feedback, can be interpreted already at 
a symmetry level, as a “closeness condition” in the symmetry space (see the 
concluding Remark). 

Our results are based on the observation that two fundamental aspects are 
involved in pattern formations: symmetries on the one side (algebraic content) 
and changes in time and space on the other side (differential content). In par¬ 
ticular, to keep account of the ‘interaction’ of both aspects, we act effectively 
with open Lie algebraic structures by using their refined structure: we need to 
introduce, a notion which generalizes the concept of a homogeneous space, i.e. 
that of an algebraic skeleton E = q(BV on a finite-dimensional vector space V, 
with g a possibly infinite dimensional Lie algebra (see Appendix 1). 
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Consider then the following infinite dimensional vector space E = g © V 

[0,0] = o, [0,0] = 0, [0,0] = 0, [0,0] =?, 

[0,0] = 0 , [- 01 , - 07 ] =0, [0,0] =0 ••• 

[ 0 , 0 ] = - 22 > 0 , [ 0 , 0 ] = 2 £> 0 , [ 0 , 0 ] = ?, [ 0 , 0 ] = - 22 > 0 , 

[ 0 , 0 ] = 0 , [ 0 , 0 ] = [ 0 , 0 ] + [ 0 , 0 ] ' • ' 

[•03, -04] = 2Z>0s, [03, 05] = ?, [03, 06] = 2 k£>0, 

[03,07] = [02,08] - [04,06], [03,08] = 0 • • • 

[04,05] = 0 , [04, 0e] = 05], [04,07] = [02,05] “ 2250, 

[04,0s] = ^>[03, 05] ' ' ' 

[05,06] = [0, [03,05]] ~ [03, [02, 05]], [05, 07] = [04, [02,05]], 

[05,08] = [04, [03, 05]] ' ' ' 

[06, 0] = 2£>0 - [01, [02,05]], [06, 0] = -2kX>0 - X>[0, [01,0s]] ' • ' 
[07, 08] = [04, [03,07]] - [03, [04, 07]] ' ‘ ' 


where D, n are real parameters and “?’ denotes undefined commutators. 

Note that the commutators [0,0], [0,0] , [03,0s] are not defined and in 
particular that even introducing new generators 0g, 0o, ip n, the algebra anyway 
does not close as a Lie algebra and is therefore an ‘open’ Lie algebra structure. 
Let us describe its refined structure. To this aim, it is also important to stress 
that many other brackets between elements of the above vector space E are not 
defined as Lie algebra brackets since they are given in terms of [0i, 0s], [02,0s], 
[03, 05]- 

The vector space V is finite dimensional and generated by [0i,0s], [02,0s], 
[03,05], 0. It has the property that each bracket of 0i, 02,03,04,06,07,08 
(freely generating an infinite dimensional Lie algebra g) with its generators is 
again in V. The Lie algebra g can be characterized as a Lie algebra of in¬ 
variant vector fields tangent to the fibers of a bundle with structure group G; 
therefore 0 can be regarded as linear differential operators by the correspon¬ 
dence 0j ipf -t— , while the Lie brackets can be interpreted as commutators 

[0i W’ = $ W^W ~ W = with the structure 

constants of g. 

In particular, we note that the commutator relations [03,07] — [0,0] = 

£*[01,05], [0,0] - 2£>07 = [0, [0,0]], [0,0] + 2 nVipg = £>[0, [0,0]], 
[0,0] + 2X>0 = [0,0] and the related [0,0] = [0, [0,0]] - [0, [0,0]] 
say that unknown commutators in the freely generated Lie algebra are related 
in such a way that their assigned relations are elements of V ; according to 
Appendix 1, the above algebraic relations define an infinitesimal skeleton. 

It is important to understand that such a feature of the commutator relations 
constitutes a skeleton’s peculiar property which generates nonlinearity. Nonlin¬ 
ear terms appearing in the ‘generated’ differential systems are dual to such rela¬ 
tions, therefore commutator relations of skeletons are the algebraic counterpart 
of nonlinear population interactions. 
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2.1 Models for pattern formation on the shells of molluscs 

In order to exploit symmetry properties, we can define an homomorphism of 
the infinite dimensional freely generated Lie algebra with a quotient Lie algebra 
by fixing unknown commutators; we can do this either by fixing the value of 
[0i,0 5 ], [02,0s], [03,05] as generated by elements in the freely generated Lie 
algebra, or by fixing the value of the unknown commutators [ 02 , 0 s], [ 03 , 07 ], 
[04,07], [06,07], [06,08]- 

Details on various choices can be found in Appendix 2. Here we point out 
the case obtained by putting 07 = 0 and 06 = 0 (since 03 ^ 0 and V 0 
this implies n = 0, i.e. the case where one of the reaction coefficient vanishes) 
we obtain an homomorphism with the Lie algebra corresponding to a group of 
Euclidean movements in the plane 

[01,02] = o, [01,05] = 0, [02,05] = 2£>02 . 

The above results seem to be in agreement with what announced concern¬ 
ing the case T> > 0 in [B], whereby, instead, integrability of reaction-diffusion 
type equations with V < 0 , i.e. diffusion constants of opposite sign, have been 
studied (such systems play a role in gauge theory of gravity [6] ; similarities with 
coupled nonlinear Schroedinger equations describing the waves propagation in 
optical fibres can be recognized, see e.g. (SI; both cases are integrable and 
homomorphisms with infinite dimensional loop Lie algebras have been deter¬ 
mined) . 

However, let us stress that alternatively we can choose 
[01,05] = A01, [02,05] = A03, [03,05] = 0 

which define an homomorphism with a finite dimensional Lie algebra with spec¬ 
tral parameter A. 

This latter case implies T> = 0, i.e. the diffusion constant of the activator is 
zero; the system would be integrable and would admit a Lax pair. In fact, the 
case T> = 0 has been related with existence of travelling waves; it appears in 
a model for pattern formation on the shells of molluscs [5)J and it is perhaps 
important to note that it is a limit case of T> « 1. In the following we shall 
consider the case 0 < V « 1 and we shall obtain the Koch & Meinhardt 
activator-substrate system from an integrability condition for a tower with a 
slightly modified skeleton. 

3 Reaction-diffusion models from integrable tow¬ 
ers with skeleton E 

We shall now explain how skeletons can ‘generate ’ nonlinear differential systems. 
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3.1 Twisted reaction-diffusion models 

For the purpose of this work, in fact, it is now a very remarkable feature that it 
is possible to obtain reaction-diffusion type models directly from the skeleton by 
using a generalization (a truncated version) of the structure equations. Indeed 
we need to introduce a way to produce (exterior) differential equations. The 
differential structure can be modelled on the assigned skeleton by its absolute 
parallelism forms (see Appendix 1); the corresponding integrability conditions 
are given by 

w fc = d£ k - , du k =rfd6 j A0* = O, (mode/); 

where 6 k are some linear differentials (horizontal 1-forms on P —>• Z) and A 
stands for the skew (also called exterior) product. 

This provides us with the differential content we are looking for. In fact, the 
forms u k can be recognized as tower forms of Cartan type ( i.e . as a pull-back 
of contact forms by a Backlund map m) if the following exterior differential 
contraints are satisfied by the 9 k : 

d6 1 = 0, d<9 4 = 0, dO 5 = 0, 

dO 2 - 2V8 2 A 9 6 - 2 V9 a A 0 7 = 0, d9 3 + 2 nV8 3 A 9 6 = 0, 
d9 6 - V8 2 A 9 3 = 0, dd 7 + V9 2 A 9 4 + 2 9 6 A 9 7 = 0, 
d8 s + V9 3 A 9 4 + 2 kV9 6 A 9 8 = 0, 

9 1 A 9 5 + V9 2 A 9 8 + V9 4 A 8 6 = 0, 9 2 Ad 5 + 8 4 A 9 7 = 0, 

9 2 A 9 8 + 9 3 A 9 7 = 0, V8 3 A 9 5 — 6» 4 A d 8 = 0, 

8 5 A 8 6 = 0, 9 5 A 8 7 = 0, 9 5 A 9 8 = 0, 9 6 A 9 8 = 0, 9 7 A 9 8 = 0. 

A solution is given by: 

9 1 = dt , 9 5 = —dt , 9 6 = ——p wdt , 9 7 = — fidt, 

V V v 

9 8 = 9 4 = ^ dx , 6> 2 = + pt^dt), d 3 = -^(udx + v x dt ), 

where /a and za are functions, depending on x and t , which must satisfy the 
following twisted reaction-diffusion equations 

pq d !dxx 2pz za 2pz = 0, v t i^xx T 2zcpiz/“ = 0 5 

here, as usual, the supscripts means partial derivatives. 

This system can be recognized as a twisted reaction-diffusion model with a 
zero basic production term , with D a normalized diffusion constant and n a nor¬ 
malized cross-reaction coefficient. It contains, as a limit feature, the particular 
case V = 0 of a model for pattern formation on the shells of molluscs |8\ 
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3.2 The activator-substrate reaction-diffusion model 

In [T5] it has been pointed out that various examples of real ecosystems, such 
as arid, wetland, savanna ecosystems, coral reefs, mussel beds, etc. can be 
described as activator-(depleted)substrate systems [3] SI In the following 
we shall obtain such type of systems by acting on the skeleton of the twisted 
reaction-diffusion system above. 

As we already explained the differential structure {i.e. the absolute par¬ 
allelism) is totally of general nature, while the somewhat ‘true’ content of a 
specific model comes from the algebraic stuctures we insert in the structure 
equations. This fact suggest the possibility of characterizing different models 
by their algebraic content. 

Let us now, in fact, consider a slight change in the algebraic skeleton, such 
as 


{ 1 P 2 , = 2nVif 3 - 22>0 2 , [V>3, i/’e] = 0, [V^, V’s] = 2nVf’ 3 . 

Such a change provides the exterior differential equation 
d9 3 - 2 kV9 2 A 9 6 - 2 kV9 4 A 9 5 = 0 , 
and therefore originates the system: 

- E>hxx ~ 2 O/i 2 - n) = 0, v t — v xx - 2k(1 - vp 2 ) = 0 . 

i. e. the activator-substrate reaction-diffusion model proposed by Koch & Mein- 
hardt (in a more general form previously also appeared in Gierer & Meinhardt 
my here again I? is a normalized diffusion constant and n a normalized cross¬ 
reaction coefficient. 

Remark 1 We stress that the latter system has been obtained by operating a 
slight change in the algebraic skeleton of the former and it is rather different 
from that one. Note, in particular, that activator-substrate systems have the 
property that both time derivatives p t and v t (activator production term and 
depletion term) are proportional to p 2 and depend linearly on v, while the 
system obtained in the subsection above is in a symmetrically twisted form. 

We see, indeed, that commutator relations of the skeleton are the algebraic 
counterpart of nonlinear population interactions. This implies that the param¬ 
eters appearing in a given model, and even the form of the model itself can be 
somewhat ‘controlled’ already at an algebraic level. This is in agreement with 
the fact that, while Lie group actions provide algebraic forms of dynamics, ‘de¬ 
formations’ of Lie algebraic structures, such as skeletons, provide the nonlinear 
content. 

Moreover, suppose a metrics could be defined on E, so that one can reason¬ 
ably think of a condition of closeness, which we write for simplicity as [^4,^7] 
— \ip 2 ,ips]- Having a glance at the skeleton structure, it is evident that such 
a condition would be equivalent to the request that D << 1, i.e. the diffusion 
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constant of the activator be much lesser than the diffusion constant of the sub¬ 
strate; therefore, we can characterize a condition for the appeareance of patterns 
by means of properties of vectors in E. 

Finally, in order to highlight and emphasize the significance of our approach, 
let us stress that the latter statement is a global symmetry-related expression of 
(and are in agreement with) empirical observations. In fact, beside the neces¬ 
sity that organisms modify their environment by inducing a long-range negative 
feedback thus allowing regular ecological pattern formations, the strenght of 
this feedback depends on the density of the organisms at large scale [15]. In¬ 
deed, approaches based only on the actual differential expression of a system, 
necessarily of a local nature, could overlook aspects related to global properties 
of the system itself such as for example, internal symmetries. 


4 Conclusions 

Our results are of general nature and in principle could be applied to other 
mathematical models proposed in various branches of biology and ecology, see 
e.g. jlO) and, for a review and further developments, El; of particular interest 
would be the possibility of application to models with delay. In fact, spatio- 
temporal pattern formation can be caused by time delay factors. The study of 
algebraic structures generating models with delay would be therefore of partic¬ 
ular interest in comparing at an intrinsic algebraic level the various approaches 
in modeling pattern formation. 

Moreover, according to the Remark above, we enhance the possibility of an 
algebraic-geometric study of the stability of the equilibrium states. In particular 
the emerging of both Hopf and Turing bifurcations depends on the parameters 
range and their mutual relations in a given model: we saw that such characters 
could be formalized already at the algebraic level in terms of the representation 
of g on V (commutator relations of skeletons). 
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Appendix 1 

Towers with skeletons 

Let us first recall a few mathematical tools constituting the background for a 
detailed treatment of which we refer to Pii and GEE; the use of the 
concept of tower with skeleton has been inspired by (and is a mathematical 
generalization of) the procedure outlined in [5], 
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An algebraic skeleton on a finite-dimensional vector space V is a triple 
(E,G,p), with G a (possibly infinite-dimensional) Lie group, E = g © V is 
a (possibly infinite-dimensional) vector space not necessarily equipped with a 
Lie algebra structure , g is the Lie algebra of G , and p is a representation of g on 
E such that it reduces to the adjoint representation of g on itself. The fact that 
E is not a direct sum of Lie algebras, but an open algebraic structure is fun¬ 
damental in order to be able to generate whole families of nonlinear differential 
systems, starting from it. 

We now consider a suitably constructed differentiable structure which is 
somewhat modelled on the skeleton above. Let us introduce a differentiable 
manifold P on which a Lie group G, with Lie algebra g, acts on the right; P 
is a principal bundle P —» Z ~ P/G. By construction, we have that Z is a 
manifold of type V, i.e. \/z £ Z , T z Z ~ V. Suppose we have a way to define a 
representation p of the Lie algebra g on T z Z ~ V, in such a way that it could 
be possible under certain conditions to find a homomorphism between the open 
infinite dimensional Lie algebra, constructed by p, and a quotient Lie algebra. 
Let us call £ the (possibly infinite dimensional) Lie algebra obtained as the 
direct sum of such a quotient Lie algebra with g. From the differentiable side, 
a tower P{Z, G) on Z with skeleton (E, G,p) is an absolute parallelism w on 
P valued in E, invariant with respect to p and reproducing elements of g from 
the fundamental vector fields induced on P. Let then t be a Lie algebra and g 
a Lie subalgebra of t. Let G be a Lie group with Lie algebra g and P(Z, G) be 
a principal fiber bundle with structure group G over a manifold Z as above. A 
Cartan connection in P of type (t, G) is a 1 form u> on P with values in t such 
that w| t p p : T p P — > t is an isomorphism Vp £ P , R*ui = Ad(g)~ 1 u for g £ G 
and reproducing elements of g from the fundamental vector fields induced on 
P. It is clear that a Cartan connection (P, Z , G, oj) of type (6, G) is a special 
case of a tower on Z. In the following, we shall be interested in the case when 
from a tower one can construct a Cartan connection by a quotienting. 


Appendix 2 

Homomorphisms with finite dimensional Lie algebras 

We can find a homomorphism with a finite dimensional Lie algebra from the in¬ 
finite dimensional open Lie algebra generated by the skeleton E given in Section 
[2] by taking the following quotient, with A being a real parameter. 

[03, Vv] = -Mi , [03, iP4 = Mr , [^ 4 , i>7\ = Ms ■ 

We get then Xtp'j = 2V^%, [ipiMe] = Mi> an d, provided that P / 0, 

bPufal = ^4 , [V>2, M] = 2P^2 + AV>3 , ['03, 05] = ^03 ■ 
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Since we also have that [04,0s] = ^[V’3> V’s] ( see the skeleton), for the consis¬ 
tency of the relations, we get in particular A = 22); thus we can write 

[07,03] = 22?0 4 , [0 3 , 0 4 ] = 22)0 7 , [04, ip 7 ] = 22?0 3 . 

It is easy to see that in this case the algebra closes as a Lie algebra 

[01, -05] = 204 , [02, 05] = 22?(02 + 03) , [03, 05] = 203 • 

Furthermore, since 07 = 08, then [06,0s] = —2k2)0s — 2)[03, [0i, 0s]] = 
[06,07] = 0 provides, by an iterated application of the Jacobi identity, the 
condition —2kD + 22? = 0. Then, if 2? ^ 0, we must have k = 1, corresponding 
to activator and substrate having the same cross-reaction coefficients. Therefore 
we get 


[02,03] = -22>0 6 , [03, 0 6 ] = 22)03 , [02,0e] = -22)02 ; 

but the consistency of [ 05 , 07 ] = [ 04 , [ 02 , 0 s]] and [ 05 , 0 s] = [ 04 , [ 03 , 0 s]] im¬ 
plies 2? = i, which corresponds to a closed Lie algebra (without a spectral 
parameter). 

Let us then consider a different closing homomorphism by setting 08 = 0 
and [03,07] = -A04, [04,07] = A0 3 . We get [03,0s] = 0, [04,06] = A0 4 , and, 
in particular, provided that 2) 0 0, 

[01,05] = ^04 , [02, 05] = 22)02 + A03 , [03, 0 5 ] = 0 . 

By the Jacobi identity, from [0 5 ,07] = —42) 2 07, and [0i,07] = 0, we get 
[04,07] = 0 which implies A03 = 0. We esclude the trivial case 03 = 0; therefore 
we let A = 0. In this case the commutation relations become [03,07] = 0, 
[03,04] = o, [04,07] = 0, [01,05] = 0, [02,05] = 22)0 2 , [03,05] = 0; the 
remaining commutators are 

[05, 07] = —42) 2 0 7 , [05, 0 6 ] = —42) 2 0 6 , [0 6 , 0 7 ] = 22?0 7 . 

[02,03] = -22)06 , [03, 0e] = 2 k2>03 , [02,0 6 ] = -22)02 , 

while we do not get a constraint on n. 
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